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1. Introduction 

Numerical simulations are one of the main tools that allow us to obtain quantitative results 
in many fields of physics, and they have played an essential part in the progress of QCD and 
condensed matter physics in recent years. There are, however, many interesting physical systems 
for which we do not have efficient numerical algorithms yet. Examples include QCD at finite 
density or with a non-vanishing 6 term. This situation has hindered progress in such fields for a 
long time, and it is thus of great interest to study novel simulation algorithms. 

In the present work we develop and test a geometric algorithm which is applicable to the two- 
dimensional antiferromagnetic Ising model with an imaginary magnetic field id (see [|T]] and [Q]) at 
6 = 71, and which solves the sign problem that this model has when using standard algorithms. 

2. The model and a geometric algorithm 

We start with the (reduced) Hamiltonian for the Ising model in an external magnetic field, 

H{{s,},F,h) = -F ^ s.Sy-^Y.'-. (2.1) 

We denote the sites of the square two-dimensional lattice by x = {ij), with i,] S {1,2, . . . ,A'^}, 
N = 2n ^fi, and the spin variables by Sx G {±1}- The sum Y,{x,y)e.'iS r^i^^ o^^i" the set 3S of all 
nearest-neighbors (x,j); F = J/{KT) is the reduced coupling between spins, and h = 2B/{kT), 
with B the external magnetic field. The total number of spins is A'^^ = 4n^, thus even, and therefore 
the quantity Q = ^Y^x^x an integer number, taking values between —N^/2 and which can 

be thought of as playing the role of a topological charge. We are interested in studying this system 
for imaginary values of the reduced magnetic field h, i.e., for h = id. 

This model suffers from a sign problem, because the weight of a configuration is not a positive 
real number, and therefore we cannot apply a standard algorithm. For the special case 6 = n we 
will show how to construct an efficient geometric algorithm that circumvents this problem. 

The partition function of the system at = TT is 

z(F, e = 7i)= £ E(.,,)e-^^.v^.v+'f ic^c 

= '^^ H "I n [cosh(F5'.^5v) + sinh(F5;t^,,)]]^52 

= E I n [cosh(F) + sinh(F)^,5,]n^z 

It is easy to see, by decomposing the lattice in two staggered sublattices, that Z(f , n) = Z{—F, n); 
therefore, at 6 = n, the ferromagnetic and antiferromagnetic models have the same partition func- 
tion, and furthermore Z is in fact a function of 

We can now expand the product inside the curly brackets, and assign a unique graph to each 
temi in the expansion: to each factor sinh(|F|)5v^v we assign the bond {x,y) £ ^ in the graph, 
and we say that such a bond is active. Every other bond is called inactive, and corresponds to 
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a factor cosh(|F|). Each subset of active bonds describes one and only one of the terms in the 
expansion. After summing over spin states most of the contributions in the expansion vanish; only 
the graphs such that every node on the lattice has an odd number of active bonds touching it give 
a non-vanishing contribution. We call such graphs admissible, and we denote by the set of all 
admissible graphs. If we consider an admissible graph, g € ^, and denote by ^b{g) the number 
of active bonds in g, by ^b{g) the number of inactive bonds in g, and by ^ the total number of 
bonds in the lattice (thus .yV = ^b{g) + ^(g)), the weight of such graph in the partition function 
is 2^sinh(|F|) '^'(^^ cosh(|F|)'^(^). Therefore we can rewrite the partition function as: 



Z{F,d = n)=2^' £ cosh(|f |)-^*(«)sinh(jf l)-^^(^) 
= 2^'cosh(|F|)''^ £ tanh(|F|)-^*(«). 



(2.3) 



The important point here is that all configurations (graphs) have positive weights, and therefore 
this representation of the partition function does not have a sign problem. Now, generalizing the 



partition function (2.3) to the case where the coupling is position-dependent, we obtain: 



Z{{F„},d = 7i)= £ <^ n ['^o^h(|F,-,|) + sinh(|F„,|)^,-^,]n^.[ • (2-4) 

{.v,=±l} l(.v,y)e,^ z ) 

From this expression we are able to build all the correlation functions for an even number of spins 
(correlation functions with an odd number of spins are automatically zero). After some calculations 
(see [^) we obtain the following expression for the correlation functions: 



C{d,F) = {s^s^,.l) 



logZ({F^,},0 = ;r) 



{F„}={n 

= ((tanh(f , (2.5) 



where ^b[g,{xi,yi}] is the number of active bonds along the straight path connecting x and x-^r 
d\. However, we stress the fact that the specific choice of the path is irrelevant, as they are all 
equivalent, as long as the endpoints are fixed. 



3. Simulation 



In order to perform calculations by means of Monte Carlo methods, we need an efficient 
algorithm to explore the space of configurations, that in the geometric representation is given by 
the set of admissible graphs, The essential ingredient is a local prescription^ that takes the 
system from an admissible configuration to another admissible configuration. The simplest change 
that one can do to a configuration is to make an inactive bond into an active one, or viceversa. 
However, applying this change to an admissible configuration will not take us to another admissible 
configuration. Let us consider instead an arbitrary square on the lattice, and consider all possible 

'Local in the sense that only a fixed number of bonds in a bounded region are changed when updating a configura- 
tion. 
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changes in the state of the bonds in the square. It is easy to convince oneself that if we start from 
an admissible configuration, the only way to anive at another admissible configuration is either 
by not changing the state of any of the bonds, or by changing all of them. Using these steps and 
the weights of the configurations derived above, it is now easy to set up a standard Metropolis 
algorithm^ . 

Our aim is to measure how the correlation functions depend on the distance d, choosing 
different antiferromagnetic couplings f < and varying the lattice volume V = L^, in order to 
determine the staggered magnetization from their asymptotic behavior. In Figs. |] and ^ we present 
some preliminary results concerning the evaluation of these correlation functions for different val- 
ues of the coupling F and various lattice sizes ranging from L = 64 to L = 1024. Simulations 
are done collecting 100k measurements for each value of F. For each run we discarded the first 
10k configurations in order to ensure thermalization. The jackknife method over bins at different 
blocking levels was used for the data analysis. 

For antiferromagnetic couplings, the behavior of the correlation functions implies that the 
staggered magnetization is nonzero, while the total magnetization vanishes, in the whole range of 
couplings that we investigated. This result is in agreement with Refs. [|l], §], and with the mean-field 
calculation of Ref. [Q]. The apparent decrease of C{d,F), starting from d ~ 10-^20, for low values 
of is probably due to the heavy-tailed probability distributions of the correlators, which is also 
the cause of their noisy behavior. In Fig. ^ we show the probability distributions of the logarithm of 
the correlators for lattice size L = 64 and F = —0.4 and F = —2.0. We can clearly notice that for a 
low coupling \F \ the values are spread in a wider range than for F = —2.0, and also that a long tail is 
developed for large distances. This makes more difficult to obtain an accurate numerical evaluation 
of the correlators. More details will be presented in [^, where also the three-dimensional version 
of the model will be investigated. 
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In order to achieve ergodicity for periodic boundary conditions, we have to implement a few, global steps (more 
details in ||3||) 
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Figure 1: Dependence of the correlation functions on the distance d for different values of the coupling F 
and for various lattice sizes L. 
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Figure 2: Dependence of the correlation functions on the distance d for different values of the coupling F 
and for various lattice sizes L. 
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Figure 3: Probability distribution of the logarithm of the correlator C{d,F) for various distances d, for a 
lattice of size L = 64, and for coupling F — —0.4 (up) and F = —2.0 (down). 
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